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The profiles of local density of states (LDOS) around single vortices and vortex pairs in type-II 
superconductors are studied taking account of the interference of quasiparticle waves experiencing 
Andreev reflection within the vortex cores, hybridization of vortex core states and normal reflection 
at the boundaries or defects. For subgap energy levels these interference effects reveal themselves 
in a nontrivial dependence of the positions of the LDOS peaks on the intervortex distance and 
sample size: the peak positions generally do not coincide with the superconducting phase singularity 
points. The LDOS proflles are calculated for three generic examples which can be realized, e.g., 
in mesoscopic superconductors: (i) vortex-vortex pair; (ii) vortex positioned near a flat boundary; 
(iii) vortex positioned in the center of a superconducting disk. The resulting evolution of the 
Andreev interference patterns could be observable by scanning tunneling spectroscopy techniques 
in mesoscopic superconductors or disordered vortex arrays. 
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I. INTRODUCTION 



Physics of vortex matter in superconductors has been a 
long standing topic of intensive research for many years. 
Over a last decade growing interest to this topic has been 
stimulated by technological advances allowing deeper in- 
sight into the structure of vortices as well as investiga- 
tion of new exotic vortex states. In particular, local 
density of states (LDOS) measurements with the help 
of scanning tunneling spectroscopy techniques have been 
proven to be an effective tool of experimental study of 
electronic structure of Abrikosov vorticesii^i^i^i^. The 
observation of the zero-bias anomaly of tunneling con- 
ductance at the center of each vortex in these experi- 
ments clearly confirmed the existence of bound vortex 
core states predicted by Caroli, de Gennes and Matri- 
con (CdGM)^. The wave functions of the subgap states 
are localized inside the vortex core because of the An- 
dreev reflection of quasiparticles at the core boundary. 
For each individual vortex the energy e(/i) of a subgap 
state varies from — Ag to -I-Aq as one changes the an- 
gular momentum jj, defined with respect to the vortex 
axis. At small energies |e| ^ Aq the spectrum is a linear 
function of /j,: e:(/x) ~ — A"^, where oj « Aq/ (kp£^), Aq is 
the superconducting gap value far from the vortex axis, 
kp is the Fermi momentum, ^ = hvp/Ao is the coher- 
ence length, vp is the Fermi velocity, and fi is half an 
odd integer. Note that hereafter wc assume the Fermi 
surface to be a cylinder and neglect the dependence of 
the quasiparticle energy on the momentum component 
along the vortex axis z considering a motion of quasipar- 
ticles only in xy plane. Introducing a cylindrical coor- 
dinate system (r, 9, z) and defining an impact parameter 
as 6 = —fi/kp = [r,ki7'] ■ Zg/kp the quasiclassical LDOS 
inside the isolated vortex core {r can be found as 

follows: N = {kp/2TTC) Sq"" 5{e ~ ujkpr sm{e - B^))dQ^. 
Here wc evaluate the LDOS summing up over the quasi- 



particle states at the trajectories characterized by the 
orientation angle dp. This expression yields a singu- 
larity of zero energy LDOS at the vortex center— i^ii^: 
N = l/(27rwrO w ^Q^jr, where iVo = (l/27r)TO/ft2 is 
a normal metal LDOS. Smearing of energy levels due to 
scattering effects leads to a reduction of LDOS peak am- 
plitude at the vortex center. However, the peak in the 
LDOS distribution survives even in "dirty" limit when a 
mean free path is smaller than a coherence length I < ^ 
(see Refs.[illl|). 

The increase in the external magnetic field results in 
the decrease in the intervortex distance which should 
be accompanied by the overlapping of wave functions 
describing the quasiparticle states bound to the neigh- 
boring vortex cores. This hybridization of vortex core 
states can perturb the CdGM spectra of isolated vor- 
tices and the circular symmetry of the LDOS peaks (see, 
e.g., Ref. [1). Note that the LDOS peak anisotropy may 
also reflect possible superconducting gap and Fermi sur- 
face anisotropy, and formation of the charge density wave 
order-ii^i^. 

The goal of the present manuscript is to show by the 
example of mesoscopic superconductor that the interfer- 
ence of quasiparticle waves trapped in neighboring vor- 
tex cores can result even in more dramatic consequences 
for the LDOS profiles around the vortex configurations. 
These consequences are not limited just to the distor- 
tion of the LDOS peaks at individual vortex positions 
caused by the concrete symmetry of the system and/or 
vortex configuration. The peak positions themselves ap- 
pear to be shifted from the superconducting phase singu- 
larity points which arc usually considered as vortex cen- 
ters. Changing the external magnetic field we can control 
the vortex configuration in the sample and switch be- 
tween different Andreev interference patterns in LDOS. 
The distinctive features of the electronic structure of vor- 
tices in mesoscopic samples with sizes of several coher- 
ence lengths are controlled by two key factors: (i) quasi- 
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particle tunneling between the neighboring vortex cores 
which becomes extremely important when the individual 
vortices merge into the giant (multiquantum) vortej*^^ 
with a winding number larger than unity M > I; (ii) 
normal scattering of quasiparticles at the sample bound- 
ary. 

The multiquantum vortex provides a simple illus- 
tration of the above statement regarding the differ- 
ence in the positions of the LDOS peaks and phase 
singularity points. Indeed, the subgap spectrum of 
such vortex contains M anomalous spectral branches 
(per spin projection)^. Each anomalous branch inter- 
sects the Fermi level and contributes to the low-energy 
DOS. At low energies the spectrum has the following 
formHOiis^. 



Ao_ 



(1) 



where the index j enumerates different spectral branches 
(1 < j < Af), ~kp^ < fij < kpi. The LDOS profile 
corresponding to the spectrum ^ consists of a set of 
axially symmetric ring structuresi^ii^. Note that for an 
even winding number the anomalous branch crossing the 
Fermi level at = (i.e. at zero impact parameter) is 
absent and, as a result, the LDOS peak at the vortex cen- 
ter disappears. The splitting of a multiquantum vortex 
into the individual vortices should lead to the transfor- 
mation of LDOS rings into a set of peaks positioned in 
the cores of individual vortices. The initial stage of this 
LDOS transformation for small intervortex distances was 
studied in Ref. [l^ within a perturbation approach. A 
nonperturbative approach which allows to describe the 
spectrum transformation accompanying the decay of the 
multiquantum vortex has been suggested in Refsj ^^i^^ . 

The second factor which is crucial in mesoscopic super- 
conducting samples is a quasiparticle scattering at the 
sample boundaries which comes into play when vortices 
approach a superconductor surface or a size of supercon- 
ducting sample is small enough. For vortices positioned 
rather close to the sample surface the effect of quasiparti- 
cle reflection at the boundary on the spectrum and total 
DOS was investigated recently in Refs. _18,, 19,2j. The 
early stage of transformation of the LDOS profiles for 
vortices approachinga flat boundary has been studied 
numerically in Ref. l2ll . 

In the present paper we analyze the effect of both the 
hybridization of vortex core states and boundary scatter- 
ing on the LDOS structures and for this purpose consider 
three generic examples: (i) vortex-vortex pair; (ii) vortex 
positioned near a flat boundary; (iii) vortex positioned in 
the center of a superconducting disk. We will show that 
the LDOS profiles can be obtained qualitatively from the 
spectra found in Refs. 17,18,19. Our qualitative consid- 
erations will be confirmed by the detailed numerical cal- 
culations. 

This paper is organized as follows. In Scc.[TT]we give an 
overview of the theoretical framework which is employed 
in this work, namely the Bogoliubov-de Gcnnes theory 



and the quasiclassical Eilenberger approach. We discuss 
the spectrum and LDOS patterns for multi vortex con- 
figurations in Sec. IIII Al and address the case of a vortex 
in a small mesoscopic cylinder in Sec. IIIIBI We give our 
conclusions in Sec. HVl 



II. MODEL AND BASIC EQUATIONS 

Our further consideration is based on two approaches: 
(i) Bogoliubov-de Gennes (BdG) equations; (ii) quasi- 
classical Eilenberger theory. The first one appears to be 
more transparent for the qualitative analysis of the inter- 
ference effects and convenient for numerical calculation of 
LDOS profiles around the vortex placed in a supercon- 
ducting disk. The BogoHubov-de Gennes (BdG) equa- 
tions for particle- {u) and hole-like {v) parts of the wave 
function have the following form: 

1 / e \ 2 

-- — p A u + Av ^ (e + ep)u, 

2ni \ c / 

^ (p+-aY v + A*u = {e-ep)v. (2) 

Here A is a gap function, A is a vector potential, p = 
^ih(d/dx,d/dy), and r = {x,y) is a radius vector in 
the plane perpendicular to the vortex axis. The LDOS 
can be expressed through the eigenfunctions of the BdG 
equation ^ in the following form (sec, e.g., Refl23): 



7V(r,e) 



\Univ)\Hie ~ Sr.) 



(3) 



where u„ (r) is electron component of quasiparticle eigen- 
function corresponding to an energy level e„ (we sum over 
both positive and negative values £„). The eigenfunction 
has to be normalized: 



(|M„(r)|2 + |«„(r)p)d2r = l 



In general case BdG equations are rather complicated. 
A simplification can be obtained if one considers a qua- 
siclassical approximation, assuming that the wavelength 
of quasiparticles is much smaller than the superconduct- 
ing coherence length (see, e.g., Ref. [231 ). Within such 
an approximation, quasiparticles move along linear tra- 
jectories, i.e. straight lines along the direction of the 
quasiparticle momentum hp = kp{cos0p,sui9p). Gen- 
erally, the quasiclassical form of the wave function can 
be constructed as follows: {u,v) = e'^^'"^{U,V), where 
{U{r),V{r)) is a slowly varying envelope function. Then 
the system ^ is reduced to the quasiclassical Andreev 
equations: 

TO V C / 

!^(ihW + -A)v + A*U^eV, (4) 

TO V C / 
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which are defined at the linear trajectories determined 
by the direction of the quasiparticle momentum \s.p. 

The quasiclassical approximation allows to develop a 
powerful method for calculation of the LDOS based on 
the solution of Eilenberger equations for quasiclassical 
propagator along the trajectories^!. For numerical treat- 
ment of these equations we follow the Refs. [8ll25l and in- 
troduce a Ricatti parametrization for the Green function. 
The essence of this method is a mathematical trick which 
allows to solve two first order Ricatti equations instead of 
fourth-order system of Eilenberger equations. Following 
Refs. HUli one can obtain: 

d 

hvF-g^a{x') + [2ojn + A*a{x')]a{x') - A = 0, 

hvF^b(x') - [2(D„ + Ab(x')]b(x') + A* = 0, (5) 
ox' 

where x' = (ki?r)/fci? = rcos{9p — 9) is the coordinate 
along trajectory, A(r) = |A|e'*, <i>(r) is a superconduct- 
ing phase, iujn — ioJn + mwpWg is a Dopplcr-shifted Mat- 
subara frequency, uon = 7rT(2n -I- 1) and 

is a gauge-invariant superfluid velocity. The LDOS may 
be expressed through the scalar coherence functions a 
and b as follows^i 

where e is the quasiparticle energy measured from Fermi 
level and is a parameter which accounts for an elas- 
tic scattering. Throughout this paper we will assume 
the simplest model^^ for the gap function distribution 
around an isolated vortex positioned at the origin A(r) = 
Ao/i(r), where 



with the core size — To describe the system of 
two singly-quantized vortices positioned at r = ±a/2 = 
±(a/2,0) we fit the gap function as follows: 

A(r) = Ao/i(r-a/2)/i(r + a/2). 

The image method allows us to reduce the problem of 
a vortex positioned at the point r = (a/2, 0) near a flat 
boundary x = to the one describing a vortex-antivortex 
pair with the antivortex situated at r = —(a/2,0). For 
the latter case we choose the gap function in the form: 

A(r) = Ao/i(r-a/2)/*(r + a/2). 

Certainly, the above simplest approximation for the 
gap function might not be enough to describe closely 
spaced vortex-vortex or vortex-antivortex pairs. How- 
ever, the resulting LDOS patterns obtained below are 



mainly controlled only by the distance between the phase 
singularity points and appear to be weakly sensitive to 
the concrete profiles of the order parameter absolute 
value. The key features of the LDOS patterns appear 
to hold even in the limit of the zero ^i, value, i.e., even 
for phase vortices without any suppression of the gap 
value inside the cores. 

Considering hereafter only the vortex configurations of 
the finite size comparable with several coherence lengths 
we neglect the vector potential which is known to give 
only a moderate renormalization of the interlevel spacing 
Lo (see Ref.i^). 

III. SPECTRUM AND LOCAL DENSITY OF 
STATES 

A. Vortex molecules 

We now proceed with the calculation of LDOS profiles 
for two systems: (i) vortex-vortex pair and (ii) vortex 
positioned near the flat superconducting boundary. The 
image method discussed in Ref. [l^ allows us to show 
that the latter case is equivalent to the vortex-antivortex 
pair in the bulk provided we choose the distance between 
the vortex and antivortex to be two times larger than 
the distance from the vortex to the flat surface in the 
original problem. To elucidate our main results we start 
from a qualitative description of the quasiparticle spectra 
following Ref. [H. 

Let us consider two vortices oriented along the z 
axis. In the plane {xy) the corresponding phase sin- 
gularity points (or vortex centers) are positioned at 
ri^2 = (±a/2,0). For quasiparticles propagating along 
the classical trajectories parallel to ki? we introduce the 
angular momenta /i = [r, kj?] • Zq = kFf'sm{9p — 9) and 
fli ~ fi — [ri,ki?] • zq defined with respect to the z-axis 
passing through the origin and with respect to the i-th 
vortex axis passing through the point , correspondingly. 
Neglecting the quasiparticle tunneling between the vor- 
tex cores the wave function can be found as a super- 
position of two states localized at different vortices and 
having close energies: £yi = —Luifi — (fci?a/2) sinflp] and 
£„2 = ~'^[a* + ikpa/2) sm.9p]. The transformation of the 
quasiclassical spectrum occurs due to the overlapping of 
the corresponding wave functions and can be analysed us- 
ing a standard quantum mechanical approach describing 
a two-level systent^^, which yields the secular equation 

{e ~ e^i){e ~ e^2) = {5ef , (8) 

and the resulting splitting of isocnergetic lines near the 
degeneracy point in the /i — plane {9p = for our case) 

UJ^l^~e± ^Lo^{kFal2f9l + [Sef . (9) 

The tunneling of quasiparticles between the vortex cores 
is determined by the exponentially small overlapping of 
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wave functions localized near the cores and results in the 
sphtting of the energy levels: 5e ~ Aoexp(— a/^). Thus, 
the estimate for the splitting ~ 5e/u} of isoenergetic 
lines in the (/x, 6'j,) -plane reads (see Eq. ©): 



(5/i(a) ~ kpS, exp — 



(10) 



In the case of a vortex-antivortex pair the non- 
interacting states localized at different vortices have the 
energies £„ = — ^^[/i + {kpa/2) sinOp] and Sav = — 
sin 6p\ . Taking into account the overlapping of 
the corresponding wave functions states we obtain the 
secular equation ([5]) with e^i = and £^,2 — Sav There- 
fore, the quasiclassical orbits near the degeneracy point 
6p — arc defined by the following equation: 



e = ±y^iLOfi)^ + (fe)2 - uj{kFa/2)0p . 



(11) 



The classically forbidden angular domain at e = has 
the width S9p — 4:Se/{u!kpa). One can assume that the 
appearance of such classically forbidden domain explains 
the deep structure in the local DOS profile observed nu- 
merically in Ref. [2ll for a vortex near the flat boundary 
of an s-wave superconductor. The classically forbidden 
angular domain results in the suppression of the overall 
DOSi^ and we show below that this mechanism is respon- 
sible for the anomalous spectrum branch disappearance 
when the vortex exits the sample. 

To study the transformation of spectrum in the entire 
range of intervortex distances we solve numerically the 
quasiclassical Andreev equations ^ to find isoenergetic 
lines in n — 9p plane. The resulting quasiclassical orbits 
for a vortex-vortex pair with the intervortex distances 
a = 2.5^ and a = 1.5^ are shown in Fig.[T^. The dotted 
lines show the orbits for a = 0, i.e., for a doubly quantized 
vortex. The splitting of the quasiclassical orbits grows 
with the decrease in the size of the vortex molecule both 
for the two-vortex and vortex-antivortex configurations. 
When the intervortex distance is of the order ^ the energy 
splitting becomes so large {Se ~ Aq) that the simplified 
equations are no longer valid and the spectrum 

is disturbed in the whole angular domain < 6p < 27r. 
One can clearly observe that with the decrease in the 
intervortex distance the CdGM orbits b = ±{a/2) sinOp 
corresponding to zero energy in isolated vortices trans- 
form into two lines = ifiQ which correspond to the 
spectrum of a doubly quantized vortex. Comparing the 
isoenergetic lines with the unperturbed CdGM ones one 
can see that the zero energy quasiclassical trajectories 
do not pass through the vortex center. It means that the 
LDOS peaks are necessary shifted from the vortex cen- 
ters as the intervortex distance becomes small enough. 
This phenomenon can be understood considering the fol- 
lowing qualitative arguments. Indeed, when the inter- 
vortex distance is large a ^ ^ the equilibrium size of 
the vortex pair is determined by the stability condition 
yielding the zero value of Lorentz force acting on each 
vortex: Fl cx: (/)o[zoiVs] = 0, where 00 is a magnetic 
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FIG. 1: Quasiclassical orbits in the ^ — 6p plane correspond- 
ing to zero energy level. Upper panel (a): two vortices at 
distances a = 2.5^ (solid line) and a — 1.5^ (dashed line). 
Dotted lines correspond to the case of a doubly quantized 
vortex, i.e. to a = 0. Lower panel (b): vortex-antivortex 
system with a — 2.5^ (solid line), a — 1.7^ (dashed line) and 
a = 1.6^ (dashdot line). 



flux piercing the vortex and Vg is a a superfluid veloc- 
ity induced in a given vortex core by other vortices and 
screening currents flowing along the boundaries. Thus, 
the local Doppler shift of the quasiparticle energy lev- 
els Ed = KkpVs appears to vanish for trajectories pass- 
ing through the vortex center. As a result, the position 
of the LDOS peaks coincide with the phase singularity 
points which are usually defined as vortex centers. How- 
ever when the intervortex distance is comparable with (or 
less than) ^ the part of the superfluid velocity induced by 
the neighboring vortex becomes strongly inhomogeneous 
inside the core region and the Vs value at a given vor- 
tex center diverges as ~ 1/a. As a consequence, the 
equilibrium size of the pair is no more determined by the 
condition of zero Vg at the vortex centers. The resulting 
nonvanishing Doppler shift for the trajectories passing 
through the vortex centers suppresses the LDOS at these 
points. Thus, for small a the LDOS peak positions do 
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not coincide with the vortex centers. 

The Fig. [TJa corresponds to the vortex-antivortex pair 
with the intervortcx distances a = 2.5^, a = 1.7^, and 
a = 1.6^. One can see that contrary to the case of the 
vortex-vortex pair the decrease in the intervortex dis- 
tance in this system leads to a rapid shrinking of the 
isocnergetic orbits. When the distance a is small enough 
all the trajectories corresponding to the zero energy are 
characterized by the ^ and 9p values close to /i = and 
Op = T^in + 1/2), where n is integer. For the intervortex 
distance a/^ = 1.5 the quasiclassical orbits completely 
shrink to these points in /i — 0p plane. In the real space 
the corresponding trajectories pass through the point in 
the middle between the vortex and antivortcx and arc 
perpendicular to the line connecting the phase singular- 
ity points. Thus, the two peaks of LDOS should shift 
from the vortex centers to the middle point and merge 
into one peak as the vortex and antivortex approach each 
other. 

For a detailed study of these effects we calculate the 
LDOS profiles for two-vortex and vortex-antivortex con- 
figurations by solving numerically the Eilengerger equa- 
tions (O and applying the expression for the LDOS 
The results are shown in Fig. [2] for the vortex-vortex 
pair with a = 1.5^ (upper panel) and a = ^ (lower panel) 
and in Fig. [3] for a vortex-antivortex pair with a/^ = 2 
(upper panel) and a = 1.5f (lower panel). In Fig. [2] 
it is shown how the LDOS profile gradually transforms 
from a two-peaked structure corresponding to two iso- 
lated vortices (upper panel) to the axially symmetric ring 
for a doubly quantized vortex. The positions of the phase 
singularity points are shown on the contour plot by red 
crosses. It can be seen that even when the distance be- 
tween the vortices exceeds the size of the vortex core the 
positions of the LDOS peaks do not coincide with the 
vortex centers. The distance between the LDOS peaks 
reduces slower than the size a of the vortex molecule and 
finally the peaks are smeared into the anisotropic ring 
structure (see Fig.l^b). The anisotropy vanishes and the 
ring becomes axially symmetric when vortices merge into 
a doubly quantized vortex. 

The evolution of the LDOS profile in a vortex- 
antivortex pair follows a different scenario (see Fig. [3]). 
When the distance between vortices is large enough there 
appear deep structures at the LDOS peaks (see Fig. [3^). 
The peaks themselves are slightly shifted towards the 
middle point at the line connecting the vortex and an- 
tivortex centers. Below the critical intervortex distance 
ac/S, ~ 1.7 the two peaks merge into one peak at the 
middle point (see Fig. [3]d). As wc discussed above this 
phenomenon is explained by the shrinkage of the quasi- 
classical orbits in the ^ — Op plane. When the vortex and 
antivortex come closer the amplitude of the LDOS peak 
is reduced and finally it disappears completely when the 
vortex and antivortex merge and annihilate. Mapping 
the solution of the vortex-antivortex problem on the one 
for a vortex near the flat surface we conclude that the 
zero bias LDOS peak appears to be positioned exactly at 



(a) 




FIG. 2: (Color online) The normalized LDOS profile in a two- 
vortex system for a/^ — 1.5 (a) and a/^ = 1 (b) calculated for 
ly = O.O6A0. Positions of vortices are marked by red crosses. 

the surface provided the distance from the vortex to the 
surface becomes less than 0.85^. 



B. Vortex in mesoscopic cylinder 

In this subsection we study another generic example 
illustrating the effect of the normal quasiparticle reflec- 
tion at the sample boundary on the LDOS profiles in the 
vortex state. We consider a singly quantized vortex sit- 
uated at the center of a superconducting disk of rather 
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(a) 




FIG. 3; (Color online) The normalized LDOS profile in a 
vortex-antivortex system for a/^ — 2 (a) and a/^ — 1.5 (b) 
calculated for v = O.O6A0. Positions of vortices are marked 
by red crosses. 

small radius R comparable with the coherence length. 
In this case a constructive interference of quasiparticle 
waves reflected from the boundary is known to result in 
mcsoscopic oscillations of energy levelsi^: 

= — cj/.t — S sm{2kpR — nfi) , (12) 

where S ^ Aq exp(-2i?/^) and |e(^)| < Aq. The re- 
markable fact is that at low energies these levels can be 
combined into two groups corresponding to odd and even 



values of the integer index ^+1/2. A spacing between the 
energy levels belonging to one group is 2uj while the en- 
ergy shift between the different groups is 2S\ cos{2kpR)\. 
Analogous shift of vortex core levels defined by the fac- 
tor exp(— i?/^) has been observed recently for a vortex- 
antivortex pair on a sphere with a finite radius R in 
Px + ipy superconductors^S,. Note that for a small su- 
perconducting disc with i? < i?c ~ (C/2) In(fcFC) 
have S ^ uj. Therefore the two groups of energy lev- 
els can be considered as continuous branches if the en- 
ergy discreteness with the scale oj is neglected. In this 
approximation one can use a quasiclassical expression 
fi = —kpr sin{9 — Op) and obtain the low energy branches 
in the form: 

£1,2 = LokpT sm{6 — Op) ± 5 cos{2kFR) ■ (13) 

Thus, for low energies e <^ Aq the quasiclassical expres- 
sion for the LDOS near the vortex center reads: 



7V(r, e) = 




The Eq. (jl4[) can be derived from the general expression 
for LDOS ([31) by setting the amplitude of the wave func- 
tion |[/„P = const which is a good approximation at the 
small distances r <^ ^ from the vortex center. Evaluating 
this expression for the LDOS we obtain: 

N{r, e = 0) = No ^ } „ , for r > rg , 

N{r,e = 0)^0, for r < ro , (15) 

where 

ro = ^ ^1 cos(2fcFi?)|e-«/« . (16) 

ukp 

Thus, the zero-bias peak of the LDOS at the vortex cen- 
ter which exists in an isolated vortex transforms into a 
ring structure which is similar to the LDOS pattern for a 
doubly quantized vortex. The only difference is that the 
radius of this ring is exponentially small comparing to the 
one in a doubly quantized vortex: tq/^ oc exp(— i?/^) <^ 
1. Therefore, this splitting of the LDOS peak in tunnel- 
ing spectroscopy experiments should be strongly affected 
by different smearing effects like disorder scattering or 
finite temperature. In particular, for the observation of 
the above effect we should assume the elastic mean free 
path to exceed the characteristic size of the LDOS ring: 
I > ^ cxp(— i?/^). Note that this condition is not very re- 
strictive because it can be fulfilled even in " dirty" super- 
conductors with I <C C because typically exp(— i?/^) ^ 1. 
The effect of smearing due to the finite temperature is 
controlled by the parameter T/S. 

To take into account the finite temperature effects 
and investigate the LDOS profiles in the whole range 
of energies and distances from the vortex core we solve 
numerically the BdG system ^ for a vortex in meso- 
scopic cylinder. We use the same method which was 
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successfully applied to study the mesoscopic oscillations 
of vortex core levelsi^, the spectrum of multiquanta vor- 
tices in mesoscopic superconductors^^ and the heat trans- 
port along vortex lines^. The resulting spectrum as a 
function of the angular momentum is shown in Fig. [4] 
which clearly demonstrates the splitting of the anoma- 
lous CdGM branch at low energies caused by the bound- 
ary effects. For large energies e ~ Aq the spectrum is 




FIG. 4; (Color online) Quasiparticle spectrum vs ^ for a 
singly quantized vortex in a disk of the radius R — 2^. The 
CdGM spectrum is shown by the solid line. We choose here 
fcfC = 200. 

not well described by the Eq. ([T2|) . 

Using the spectrum and the expression ([3]) for the 
LDOS we can find the local differential conductance as a 
function of voltage V: 

dl__(dl\ f+°° N{r,e) df{e-eV) 
dV^\dv)^J_^ No dV ^' 

where {dI/dV)N is a conductance of the normal metal 
junction and f{E) = [1 + cxp{E/T)]~^ is a Fermi dis- 
tribution function. The typical plots of the differential 
conductance are shown in Fig. [5] for different disk radii. 
The zero bias peak at the vortex center is clearly seen 
for a large disk radius when the normal scattering at 
the surface can be neglected. The decrease in the disk 
radius results in the shift of this conductance peak to 
higher voltages and formation of ring structure of local 
zero-bias conductance with finite radius r ~ rg. 

Note, that normal scattering of quasiparticles at the 
boundaries can be also important for multiquanta vor- 
tex configurations, e.g., vortex pairs, trapped in rather 
small samples. In our consideration of the LDOS profiles 
for a vortex-vortex pair in the previous subsection we 
have neglected the normal scattering at the sample sur- 
face assuming the vortex centers to be situated rather far 
from the boundaries. However, from our consideration of 
the vortex - antivortex problem one can expect that the 




FIG. 5: (Color online) Three-dimensional plot of the nor- 
malized local differential conductance dl /dV as a function of 
voltage {eV) and distance from the vortex center (r) in disks 
with a large radius R = 7^ (a) and small radius i? = 2^ (b) 
for T/Ao = 0.005. We choose here kp^ = 200. 



effect of boundary scattering on the LDOS becomes im- 
portant only when the distance from the vortex center 
to the sample edge is comparable with ^. In this case 
the LDOS profiles in a vortex pair are affected by both 
the hybridization of the CdGM states and normal scat- 
tering. Considering the particular case of a doubly quan- 
tized vortex one can expect the boundary scattering to 
result in the spectrum oscillations (see Fig. 6) similar to 
the ones in a singly quantized vortex. These oscillations 




-200 200 1^ 



FIG. 6: (Color online) Quasiparticle spectrum vs /i for a dou- 
bly quantized vortex in a disk of the small radius R = 2^. 
The spectrum for a doubly quantized vortex in a disk of the 
large radius R — 7(, is shown by the solid line. We choose 
here kp^ ^ 200. 
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should cause the broadening of the ring in the LDOS 
pattern with the typical broadening scale rg. The same 
effect should result in the broadening of the anisotropic 
ring structure for a vortex pair of a finite size. 

IV. CONCLUSION 

To summarize, we have investigated the distinctive fea- 
tures of the LDOS profiles specific both for the vortex 
states in mcsoscopic superconducting samples and for 
vortex arrays in bulk superconductors. These profiles 
are shown to be strongly affected by the Andreev in- 
terference effects and the normal scattering at the sam- 
ple boundaries. In particular, we find that in contrast 
to the case of isolated vortices far from the boundaries 
the positions of peaks in the LDOS profiles in vortex 
configurations in small size samples are shifted from the 
phase singularity points. Taking a generic example of the 
two vortex system, we have considered the evolution of 
the LDOS profile which accompanies the merging of two 
vortices and appearance of the doubly quantized vortex. 
In this case the distance between LDOS peaks reduces 
slower than the intervortex distance a. Finally, when a 
is smaller than the vortex core size the two-peaked struc- 
ture of LDOS transforms into an anisotropic ring. The 
anisotropy vanishes when the vortex positions coincide 
and a doubly quantized vortex is formed. 

Earlier (see Ref. [Tsl ) it was shown that the spectrum 
of a single vortex placed near the flat surface is trans- 
formed analogously to the vortcx-antivortex system and 
when the distance is of the order of the vortex core 
size, the interlevel spacing in the vortex spectrum be- 
comes larger than the CdGM value. This effect was 
argued to lead to the disappearance of the anomalous 
spectrum branch when the vortex approaches the sur- 
face. In the present paper we have confirmed this pre- 
diction by calculating LDOS profile for the entire range 
of distances from the vortex to surface. We have found 
an amazing effect for such system: as the vortex ap- 
proaches the surface the zero-bias peak of LDOS shifts 
from the vortex center to the surface point positioned 
at the minimal distance to the vortex center. For the 
particular vortex core model that we have used in this 
paper the shift of the LDOS peak occurs when the dis- 



tance from vortex to surface is 0.85^. Experimentally 
the LDOS peaks positioned very close to the supercon- 
ductor boundary have been recently observed in vortex 
state of tungsten (W) based thin films by scanning tun- 
neling microscopy/spectroscopy techniques^. Note that 
in the latter experiment LDOS vortex peaks have been 
studied near the superconductor/normal metal (Au) in- 
terface. We suppose that the normal scattering of quasi- 
particles at such boundary caused by either the Fermi 
velocity mismatch or some surface barrier can result in 
the transformation of the LDOS profiles analogous to 
the one discussed above for a superconductor/insulator 
boundary. 

Also we have investigated a measurable consequence of 
the giant mcsoscopic oscillations of vortex core levels in 
finite size superconducting systemsi^. Considering the 
simplest case of a vortex positioned at the center of a 
superconducting cylinder with a radius R we have shown 
that the zero-bias peak of the local differential conduc- 
tance at the vortex center transforms into a ring structure 
of the radius of the order of ro ~ ^e~'^^^^. We expect 
that the unusual behavior of the LDOS profiles which we 
have investigated could be observable in scanning tunnel- 
ing microscopy/spectroscopy experiments in mcsoscopic 
superconductors which are now in the focus of interest 
(see, e.g., Refs. [s^P^'s^). 

Certainly, the Andreev interference patterns discussed 
above should be smeared by the disorder effects which are 
controlled by v parameter within the Eilenberger theory. 
The resulting smearing of the LDOS profiles can be ex- 
tremely important, e.g., for amorphous superconducting 
samples with rather small mean free path. Still we expect 
that the shift of the LDOS peaks from vortex singularity 
points caused by the Doppler effect will be observable 
even in this limit of strong disorder. 
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